function Fig1_b()

data_name = '22';
path = 'raw data/';
base_name = 'GBsweep_';

data_c = '';
data = '';
menu = '';
data_ind = 1;
tmp_ind  = strsplit(data_name,',');

for i = 1:length(tmp_ind)
    list{i} = [base_name , num2str(str2double(tmp_ind{i}),'%.3d') ];
end
for i = 1:length(list)
    for ind = 1:10000
        try
            [data_c{data_ind}, ~] = load_3d([list{i} '_' num2str(ind,'%.4d') '.3dc'], path);
            data_c{data_ind}.gate_slow = data_c{data_ind}.gate_slow(ind);
            menu{data_ind} = data_c{data_ind}.gate_slow;
            data_ind = data_ind + 1;
        catch
            break;
        end
    end
end
[~,I] = sort(cell2mat(menu));
data_c = data_c(I);
data.data_didv = data_c;
data.ind = length(menu);
data.gate_slow = cell2mat(menu(I));
n_gate_slow = length(data.gate_slow);



for i_m = 1:1:n_gate_slow
    tmp_didv = data.data_didv{i_m};
    n_vb{i_m} = length(tmp_didv.gate_fast);
    vb{i_m} = tmp_didv.gate_fast;
    n_bias = length(tmp_didv.bias);
    bias = tmp_didv.bias;
    ref = 1:n_bias;
    for i_b = 1:1:n_vb{i_m}
        ind = ref(abs(bias)<1e-4);
        didv = tmp_didv.map(:,i_b);
        didv = movmean(didv,5);
        map_zbdidv(i_m,i_b) = didv(ind) *25; %2mV oscillation, 50mV sesnsitivity
    end

end


Fcolor='k';
Bcolor='w';
Fsize=18;
Bsize=14;
f = figure();
f.Color=Bcolor;



imagesc(vb{1},data.gate_slow,map_zbdidv(1:1:end,1:1:end));
ax = gca;
ax.YDir='normal';
ax.XColor=Fcolor;
ax.YColor=Fcolor;
ax.FontSize=Fsize;
xl=xlabel('Vg (V)');
xl.FontSize=Fsize;
yl=ylabel('Vm (V)');
yl.FontSize=Fsize;




cMap=flipud(getColor('Red'));
cMap=cMap(1:170,:);

cMap = flipud(cMap);

colormap(cMap);
ax.CLim=[0.1 0.5];

end